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Abstract. As can be inferred from present experiments in ultracold gases, the 
scattering length is a quantity that determines the thermodynamic state of the gas. As 
such, there exists a conjugate thermodynamic to it. Here, we show that the recently 
introduced contact is the conjugate of the inverse of the scattering length. We find that 
this identification allows for a derivation of essentially all the known results regarding 
the contact. Using the mean-field theory for the Bose-Einstein (BEC) to Bardeen- 
Cooper-Schriefer (BCS) crossover, we also find that the contact is proportional to 
the square of the gap. We analyze in detail both a homogenous balanced mixture of 
fermions and its inhomogenous counterpart in a harmonic trap. 

PACS numbers: 67.85.Lm Degenerate Fermi gases, 67.85.-d Ultracold gases, trapped 
gases, 03.75.Hh Static properties of condensates; thermodynamical, statistical, and 
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1. Introduction 

The experimental realization of the crossover between a molecular Bose-Einstein 
condensate (BEC) and an atomic Bardeen-Cooper-Schriefer state (BCS) in fermionic 
alkali vapours, such as ^Li and ^°K, is one of the most fundamental and far reaching 
results of the current ultracold quantum fluids research. The route to this crossover of 
macroscopic states started with the early observations of molecular BEC[T1 |2l [Sj HI |5] 
continuing to the observation of actual Fermi pairing in the BCS side|6l [TJ El E] that 
directly verified the crossover, and in the meantime, the superfiuid nature of the quantum 
fluids was established fTTl [T2| [T3] . Furthermore, a whole class of spectroscopic 
techniques have been developed to enquire about the elementary excitations spectra [T^ 
[T5l [T6t [T7] which not only should verify the superfiuid characteristics of the state but 
should also shed light on their true microscopic details. Of additional relevance, and the 
matter of the present article, is the study of thermodynamic properties [T8 | [T9 | [20 t |2 H [22] 
and, in particular, of direct measurement of the contact by Stewart et al.[23], a 
thermodynamic variable we discuss below. We refer to the review by Giorgini et al.|24] 
for further discussion and additional references. The BEC-BCS crossover occurs as a 
consequence of the change of the atomic scattering length by means of external magnetic 
fields, with the concomitant presence of a two-channel Feschbach resonance |251 126] . 
that permits a switch between the two extreme states. This phenomenon, as originally 
predicted by Eagles [2^ and Leggett[28]. can also be theoretically studied considering 
a one-channel potential resonance that allows for the change of sign of the scattering 
length, indicating the existence of the molecular BEC and atomic BCS. Among the 
many studies sparked by these experiments, we are concerned here with the appearance 
of a new variable called the contact. This was originally introduced in a series of papers 
by Tan[29l [30l [31], and further analyzed by others [321 [331 EH [35], and it has proved to 
be of relevance in the determination and understanding of BEC-BCS crossover. 

The contact, as introduced by Tan [29 j. tells us about the asymptotic behavior of 
the wavevector occupation number of fermions n^, and it is defined as 

C = lim fcVfc. (1) 

That is, since the contact approximation of the interatomic potential u{r), implemented 
by replacing 

, . Anh^a^,^ 

u(f) ^ 5(f), (2) 

m 

with a the scattering length, necessarily introduces difficulties as f — )■ 0, Tan found 
that this produced a divergent but well characterized behavior of the occupation 
number as — ?■ oo. This, in turn, opens a door to "correct" the well-kown ultraviolet 
divergences produced by the contact approximation ([2|). Tan called the limit ([1]) the 
contact. Further analysis [29] showed that an adiabatic change of the scattering length 
a, relates the contact to the energy E of the system by means of the relationship. 




C=-^y{^] , (3) 
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where V is the volume of the system, and the derivative is taken at entropy 5* constant. 
Here and throughout, we shall use r/ = 1/a as the inverse of the scattering length. 

In this article we point out a further and very simple connection between the 
scattering length a and the contact C, a relationship that seems to have gone unnoticed. 
We shall argue that since a, or better 77 = 1/a, is a bona fide intensive thermodynamic 
state variable, it must have an extensive conjugate variable. We show here that this 
variable is the contact C (times the volume V). We shall see that expressions ([1]) 
and ([3]), among other properties, follow from an analysis of general thermodynamic 
considerations and of this identification. To illustrate our results, we shall analyze 
in detail all the thermodynamic predictions of the mean-field theory [28] of the BCS- 
BEC crossover at zero temperature, T = 0. As we shall see, the results drawn 
from this approximate theory are in reasonable agreement with both other more 
precise theories|Mj and recent experiments [9l [23]. This work, we believe, adds to the 
comprehension and potential use of the contact variable. 

2. The thermodynamic origin of the contact 

The observation we made here is based on the fact that any mechanical variable of a 
system that can be externally varied, and that can change the current thermodynamic 
state of the system, is in itself a thermodynamic state variable. As an example, any 
adiabatic change of an internal parameter of the Hamiltonian of the system (keeping, 
say, the volume and number of particles constant) necessarily results in an increase or 
decrease of the temperature of the system; it is thus a thermodynamic parameter that 
determines the macroscopic state of the system. Let us therefore take = 1/a as a 
thermodynamic variable. Using the inverse of a is a matter of convenience. By its mere 
nature, r] is an intensive variable, namely it does not scale with the size of the system. 

Since most of the theoretical analyses of many-body systems is performed in the 
grand canonical ensemble, let us first consider the grand potential f2 of a homogenous 
gas. Such a function is a fundamental relationship that yields all the thermodynamics 
of the system [36] in terms of the following variables fl = fl{V,T, n,ri), where V is the 
volume of the system, T its temperature, fi the chemical potential, and we have assumed 
an additional dependence on rj. Being a function of (V, T, fi, rj) this is equivalent to state 
that 



which amounts to identify (or define) the entropy 5*, the pressure p, the (average) number 
of particles A^, and an additional variable, the contact C, as partial derivatives of fl. In 
particular, the contact is 



Thus, just as {p, V), (T, S), and (/i, A^), {rj, C) are a pair of conjugate variables, rj being 
intensive implies that C = C{V,T, fi,ri) is an extensive variable. As shown below, it is 



dn 



pdV - SdT - Ndfi - Cdr] 



(4) 




(5) 
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related related to Tan's definition of the contact C as C = {h'^/4:7im)CV. That is, Tan's 
C is an intensive variable, a contact density (or contact "intensity"), and in this case is 
a function C = C(T,fi,ri). 

The relationship to statistical physics serves to find further properties of C and 77. 
Recall that in the grand canonical ensemble, the state is given by the density matrix 

where is the number of particles operator and H the Hamiltonian of the system. This 
Hamiltonian depends parametrically on the scattering length a, as we further specify 
below. S is the grand partition function, 

S = Tr exp(-^iV - j^H). (7) 

The grand potential is i7 = — /cTlnS and the entropy is given by 

S = — /^^Tr p Inp. (8) 
From this we obtain, by substituting ([6]) and ([7]), 

TS = E - i^N -n. (9) 

where E = (H) is the average energy. The expression for S allows to solve for 
E = Q + TS - fiN, that with the use of yields, 

dE = TdS - pdV + fidN - Cdr]. (10) 

From this expression one finds that the contact is also given by. 



^ U^/5y,iv' ^^^^ 
which is the known expression ([3]) also derived by Tan|29]. It is the so-called adiabatic 
definition of the contact. It is clear that if one considers the Helmholtz free energy, 
F = E — TS, the contact is now C = —{dF/dr])Ny^T, and so on for any other 
thermodynamic potential. 

We now specialize to a fermionic many-body system with two-body interactions, 
and assume that the collisions can occur only by one open channel. Below we comment 
on the case of two-channel collisions |33 1 IM] . In general, any real interatomic potential 
u{r) depends at least on two parameters, say the depth of the potential Uq and its range 
tq. At low energies or temperatures, both parameters can be replaced by the scattering 
length a, and accordingly, the potential is replaced by the so-called contact potential, 
as given by eq.(l2]). As it is well known this approximation is plagued with ultraviolet 
divergences and several schemes have been developed to repair them. Incidentally, 
Tan's work on the contact variable is another instance to deal with this problem. As 
we shall return below, the replacement of the several parameters of the interatomic 
potential by solely the scattering length has very strong implications on the form of the 
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thermodynamic potentials. In any case, the Hamiltonian in this approximation may be 
written as 

H=-T. ^i^^'^^ + (12) 

where we have aheady considered a mixture of fermions, say, in different internal states 
cr = 1 and 2, ifji and ip2 being their field operators. The gas is assumed to be confined 
in a rigid box of volume V and it is, therefore, homogenous. Using the results above, 
eqs. ([5]) and ([7]) with Hamiltonian ( !T2|) . one immediately finds that the contact variable 
is, 

C = ^ / (aViV2V'2V'i)^^'^- (13) 
m J 

where the average can be taken in any ensemble, even though we obtained it in the 
grand ensemble. This result was originally derived by Braaten and Platter [32|. These 
authors argue that a way of avoiding the ultraviolet divergence in f|T3l) is by replacing a 
by g{a) = a/(l — 2aA/7r) with A a cutoff parameter. Moreover, after a careful analysis, 
they find that the above value of C is precisely the original Tan definition of the contact 
C = C/V as given in equation ([1]), namely, as the coefficient of the /c~^ tail of the 
number distribution. 

If the gas is confined by an inhomogenous potential, as in the current experiments 
with ultracold gases, say in a harmonic trap [37], 

Kxi(r) = ^rnu^r^ (14) 

one must add to the Hamiltonian ( JT2l) a term Kxt(^ So- V'i^o-- This addition is by 
no means trivial since it modifies the thermodynamic description of the system. As 
discussed in Refs. [38| l39t l40| |41] |42] . among others, the thermodynamic limit is now 

— )■ oo — )■ 0, with Nu}^= constant. This yields the "generalized" volume V = 
as an extensive thermodynamic variable that replaces the usual volume of a rigid-box 
potential. In addition there appears a conjugate "generalized" pressure V that replaces 
the hydrostatic pressure p. The details of these variables as well as their physical 
interpretation has been given at length in Refs. |10l HU US]. Here, suffice to say that 
the grand potential of the trapped gas is a function Q = Q{V,T, fi,ri). The contact is 
simply, 

^ (^^)v,T,,. ^^^^ 

and it is now a function of (y,T,fi,ri). Similarly to the homogenous case, eq.Q, one 
has the Euler relationship, Q = E ~ TS — ^N, and the generalized pressure may be 
identified as i7 = —VV, namely as the grand potential per unit of generalized volume. 
We shall use this result below. 

With the use of the local density approximation (LDA) one can calculate all 
thermodynamic variables in the trapped gas. In particular, following the procedure 
of Ref. ^2], one can show that if one has knowledge of the contact density c/i(T, /i, 77) = 
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C(y, T, fi, ri)/V of the homogenous case, the contact in the trapped gas can be obtained 
by LDA prescription, 

C(V,T,/i,r7) = I c^(T,/i- V;,i(f),r7)dV. (16) 

This result will be used in section 4. 

If one considers two-channel collisions, such as those involved by a Feschbach 
resonance, Zhang and Leggett[33] and Werner et al.[34j have shown that the number 
of bosonic molecules in the closed channel is also related to the contact C, up to a 
quantity that depends crucially on the two-body interactions involved in the resonance. 
Werner et al.[3l] find very good agreement with the experimental results by Partridge 
et al. [9]. Although we shall not consider this case here, it is of interest to mention that 
a calculation of the contact with a single channel, as performed in the present article, 
gives reliable information on a quantity that also depends on the details of two-channel 
collisions. Below we return to this point where we discuss the "universality" of the 
obtained thermodynamics. 



3. Reduced variables and "universality" 

So far, we have just shown how all the known results concerning the contact follow 
directly from the identification of = 1/a and C as conjugate thermodynamic variables. 
Here, we would like to make some general considerations regarding the functional 
dependence of the different variables by assuming the contact approximation of the 
interatomic potential, eq.([2]). That is, we assume that the interatomic interaction 
depends solely on the scattering length a. 

Let us analyze first the homogenous case. In the grand canonical ensemble 
the independent variables are (V, T, /i, t]) with only two parameters, Planck constant 
h and the atoms mass m. We note that the number of particles is a dependent 
variable = N{V,T, fi,ri). Thermodynamics ensures that one can single- value invert 
II = fi{n,T,ri) with n = N/V the homogenous particle density. Thus, any other 
thermodynamic variable can be written in terms of (A^, V, T, t]), for instance, the contact 
C = C{N, V, T, rf) or the energy E = E{N, V, T, rj). Now, as long as the particle density 
is never taken as zero, n is always positive. Therefore, one can use h, m and n to 
adimensionalize or reduce all the other variables. Given the fact that we are dealing 
with fermions, it proves convenient to use the Fermi energy and momentum 

ep = ^ and kp = (37:^71)'/ \ (17) 
2m 

and the mass m, as the quantities to define dimensionless reduced variables. We shall 
place a "tilde" ~ on top to denote reduced variables in the homogenous case. For 
instance, for intensive variables, 

u ~ kT rt . . 
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For extensive variables we will take them as per particle, namely, 
n E Ckp 

Nep Nep Nep 

and so on. We note that the number density becomes a number. 

Hence, since extensive quantities are of the form X = Nx{n,T,ri) and intensive 
ones ^ = ^{n,T,r]), then, it must be true that all the reduced variables depend only on 
T and rj, namely, /i = fi{T,ri), e = e{T,fi), c = c(T,fi), and so on. This indicates that 
the dependence is "universal", as in a law of corresponding thermodynamic states [36]. 
However, we want to make this point clear: this dependence is only possible because the 
Hamiltonian has only two parameters, h and m (a is a variable not a parameter). If the 
interatomic potential had one more parameter in addition to a, then this "universal" 
dependence will not be possible. That is, an explicit dependence on the density will 
still occur. Recently, Ho|l3] introduced a Universality Hypothesis in which, near the 
unitarity limit \a\ — )■ oo or 77 — )■ 0, all interatomic potential parameters become irrelevant 
except the scattering length a. This implies, as we see here, that all thermodynamic 
variables, for fixed values of T and fj, behave, up to a numerical constant, as if the 
system were an ideal Fermi gas. Here, by assuming at the outset that the only relevant 
interatomic potential quantity is the scattering length, we find that such a hypothesis 
is always valid. In real systems, and away from unitarity, one expects a dependence, 
maybe weak, on other interatomic potential parameters. In this connection is important 
to mention the recent work by Zhang and Leggett[33]. where they argue that near 
unitarity, where \a\ is very large and becomes the dominant length, the many-body 
physics is dictated by the contact potential and, thus, its thermodynamic contribution 
is actually universal, thus lending validity to Ho's hypothesis. Zhang and Leggett, 
though, are very explicit in the fact that certain quantities are also influenced by the 
details of the two-body physics, such as the number of bound pairs or molecules in the 
closed channel of a Feschbach resonance. As it turns out, the number of bound pairs in 
the closed channel is found to be directly proportional, both, to a quantity that depends 
on the details of the two-body interatomic potentials and to the contact C, that being a 
many-body thermodynamic quantity is independent of the details of the collision. Thus, 
two-body physics matters, but a one-channel calculation of thermodynamic properties, 
such as the contact, suffices to determine the many-body contributions. 

For the case of a gas trapped in a harmonic potential the situation is very similar. 
First, we note that the number of particles is = N{V,T, fi,7]) with V = 
Again, we solve for the chemical potential fi = fi{N/V,T,ri), and express all other 
thermodynamic quantities as functions of (A^, V, T, 77), say C = C{N,V,T,ri). Similarly, 
we introduce the Fermi energy and momentum in the trap as. 

Bp = hiSNtu'Y^' and Kp = [^^^ ■ (21) 
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We see that the combination Nu^ = N/V is a generahzed particle density for the 
harmonic trap. We define dimensionless reduced variables in the trap, denoted with a 
"hat" A, in an analogous way to the homogeneous case, 



t^ rf. kT ^ ri 
12 = — T = — r]= — 

£p £p Kp 

for intensive variables, while for extensive ones 

n ^ E ^ Ckp 

CO = e = c = , 

Nep Nep Nep 

and so on. Again, any "hat" variable thus defined depends only on T and fj, i.e. 
f = f{f,v). 



(22) 



(23) 



4. Mean-field BEC-BCS crossover at T = 

Further explicit results can now be obtained by appealing to the mean-field BCS-BEC 
theory [271 [28] at zero temperature T = 0, as clearly described by Leggett[28]. Within 
this approximation, one finds that the grand potential Q is given by. 



(efc - yu) - J{ek - /i)2 + A 



k 



m 



(24) 



where = h^k'^/2'm and the 
transcendental equation. 



gap 



2 Ve, AtiK' 
A is a function of (T, /i, 77) given by the 



-E 



(e,-/i)2 + A2 



1 



m 



2Tch 



(25) 



These two equations yield, as expected, Q = Q(y,^,ri). The thermodynamic limit is 
assumed to be used in the form J2k ~^ ^ I (27r)'^ / S'k. Few words are in order here. This 
form of the grand potential is obtained by considering that the ground state of the two- 
species fermion mixture is given by the BCS variational state[l9]. And most importantly, 
we recall that we have already included in ( 12^ . by hand, the "counterterms" needed 
to avoid the ultraviolet divergences that the contact approximation ([2]) necessarily 
introduces. Strictly speaking, the mean-field approximation is only valid in the weakly- 
interacting regime, which requires that N\a^ jV ^ 1. That is, it should only be valid in 
the limit of very large \'!]\] this immediately prohibits to use the theory in the strongly 
interacting crossover region 77 ^ 0. Nevert her less, as we shall see below, one finds 
results that agree reasonably well with current experimental results [HI [23] • We have not 
considered finite temperatures here since, on the one hand our purpose is to exemply the 
main results, and on the other, it is well known that this approach fails in considering 
correctly all the thermal excitations |44[ [351 [l6l [ITJ [18] . 

We are now in a position to calculate all the thermodynamic properties of the 
system by using as given by ([24]) and f[25|) . as a function of (V, /x, r/) and by recalling 
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that S, p, N and C are given by its partial derivatives. We are particularly interested 
in N and C. Recall that the number of fermions in each species is N/2. One finds, 



1 _ - 



(26) 



and 



C ^ V£^A^ (27) 

This second expression is an important result of this article, since one finds that the 
contact per volume C/V is, apart from constant factors, the square of the gap A. From 
the expression for N, eq. fl26|) . one can infer the occupation number of fermions (of any 
species) in state k, i.e. Uk- One sees that in the limit /c — )■ oo, the quantity k'^rik converges 
to 4:7TmC/hV = C as consistently given by (l27Il and showing, within this model, the 
equality predicted by Tan, see eq.([I]). Because of this identity, measurement of the 
contact [9l [23] may lead to an alternative determination of properties of the excitations 
energy spectra of these systems. 

Before studying some properties of the contact, from the perspective of eq.( l27Il . we 
would like to point out that if we returned to the many-body description of the interact- 
ing fermion gas, in terms of the creation and annihilation operators of particles in states 
k, and then implemented the contact approximation of the interatomic potential, eq. 
([2]), one would find that the sums are plagued by ultraviolet divergences. As Tan has 
suggested, a way to regularize the theory would be to introduce a kind of counterterm 
in the form (n^ — C/k'^). We find that in the mean- field BEC-BCS theory one can also 
regularize the theory at that level, namely, before implementing the variational scheme, 
and the counterterm is precisely Tan's prescription. At least at the level of a mean-field 
theory, we consider that this closes the circle of identifying the contact variable C as the 
thermodynamic conjugate of the inverse of the scattering length t] = 1/a. 



4.1. Thermodynamics of the homogenous gas 



We first rewrite the three main equations, f l25|) . f l26|) y fl27j) . in terms of reduced variables, 
see eqs. (|T8l) and (|T9l). 



V 



1 



TT JO 



X 



'/^dx 



and 



4 Jo 



37r 



[X - /i)2 + A2 ^ 
X — fl 

[x - fiY + A2 



(28) 
(29) 
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Figure 1. Reduced contact c vs fj and chemical potential /2 vs 77, for a homogenous 
fluid, r) = 1/kpa. 

The solution to equations ( l28l) and ( l29l) yield ft = jiirj) and A = A (77), as expected, 
and therefore, c = c{ri) as well. Fig. [1] show fi vs fj and c vs fj. At unitarity 
fj = one finds fl{0) ~ 0.59 and c(0) ~ 0.54. Typically, the value of the chemical 
potential at unitarity is related to the universal /3 parameter [50] as /i(0) = 1 + (3. From 
experiments [51 ttH [191 |20l [22] and quantum Monte Carlo simulations [SH [521 [531 [51 [55] 
as well as field theoretical calculations [561 [HI], the accepted values are (3 ~ —0.64 to 
-0.49 and A ^ 0.44 to 0.54. The mean field values are (3 ^ -0.41 and A ^ 0.68. It is 
well known that mean field is a bit off at unitarity, but as we see, such a simple theory 
allows for a full determination of the whole crossover yielding even semi-quantitative 
agreement with experimental data. 
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Figure 2. Reduced pressure p vs f] for a homogenous fluid, fj = l/kpa. We have 
plotted Stt^p to make direct comparison with the reduced pressure in the trap, see 
Fig El See text for details. 



One can also work out the asymptotic results at BCS, rj — )■ — oo, and at molecular 
BEG, fj — 7- +00. These are, BCS, 

/i ^ 1 , c^ 247re"^ e''^ as fj -oo, (31) 

and molecular BEC, 

-221^^ 

a Ri —ri H , c~2r7 as r7-> +oo. (32) 

37r rj 

Restoring the units of the first equation, it yields to lowest order at BEC, /i ~ — ^^/2ma^, 
half the binding energy of one molecule [28]. 

To complete the thermodynamics we can calculate the pressure p = —Q/V. 
However, to make comparisons with the trapped case, we shall instead calculate the 
reduced grand potential u = —Sn'^p. This can be found directly from eq.(l24ll in reduced 
units, 

3 



u = - I x^^'^dx 
4 Jo 



(x-/i)-J(x-/i)2 + A£-^c 



+ dfi. (33) 



Since ft and c are functions of fj, so is u. The result is shown in Fig. |2l Although 
not shown here, the only remaining quantity is the internal energy. This is given by 
e = u + p.. 

4-2. Thermodynamics of a gas in a harmonic trap 

This case can be found resorting to LDA. First, we restore the units of the homogenous 
case to find the equation for the density profile in the trapped case, that is, we write 
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fi = fi{n,ri) = epjl{ri/kp) and implement LDA: Take fi ^ fi — mcj^r^/2 and n — )■ p{r), 
the latter being the density profile. One finds, 

Then, the number of particles in the trap is found by integrating the density profile, 

N = J p{f)d^r. (35) 

For the calculation of the contact, one must be careful to realize that one needs first the 
homogenous contact density, namely, Ch = C/V. This is, 

c.(n,,) = ngg(^^-^). (36) 

Then one implements LDA to find the contact in the trap, 

C = j cM^,ri)dh. (37) 

Solution to these three equations yield p = p{f; p,ri,V), N = N(y,p,ri) and 
C = C(y,p,ri). To actually solve those equations one assumes N/V = Nu^ fixed and 
solves separately the cases rj <0 and rj > 0. This last separation is useful since one can 
show, using the form of p as given in Fig. [H that for fixed < 0, as /i — )• 0"*", n — )■ and 
remains zero as p becomes negative, while for fixed r7>0, n— j-Oas/x— )■ —h'^rf /2m. 
This is shown in Fig. |3l This separation allows to determine the range of values of the 
density profiles in the trap. That is, one can see that in general p{r) ^ as r — )■ r^ax, 
where Vmax is given by, 

1/2 

r —I J '/ — " ,„„x 

I max — \ f ^ ^2 o\l/2 _ • y^^j 



if ry<0 



Using now reduced "hat" variables for the trap, see eqs.( 12^ and (1251) . the equations 
to solve, flMj) and (|35|) . become, for the density profile. 



^ r■-^^2,r.^^l/3 ' (39) 



\ /i - (/i + r)2) x2 j V" ^-^-^n V(37r2p(x))i/3^ 
for the total number of atoms, 

l = 967r<^ \ \ x'^ p{x]p,fi) dx, (40) 

and for the contact, 

32 f /i3/2 ] A , ^ .4/3 ^ / fj 



c 

TT 



ip + ef' I io ('"'^("^) ' i (3.2p(x))V3 j • (41) 

In the above three equations, the upper term in braces is for 17 < and the lower one 
for fj > 0. Note that the profile becomes p = p{x;fi), with < x < 1. It is also clear 
that fj is the only independent variable, namely, p = p{fi) and c = c{fj). It is interesting 
to point out that the previous equations (!39|) - fHT]) are valid in general, not only at 
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mean-field level, since their derivation requires LDA and the boundary condition (l38l) 
only. Mean-field enters as the model for fi{fj) and c{fj). 

Equations fl39|) - PT|) can be solved numerically, with the input from the 
homogeneous solution fl{fj) and c{fj), and their solution /i vs 17 and c vs r) is shown 
in Figs, m One can find special and asymptotic behaviors. At unitarity fj = fj = 0, and 
one finds 

m = (/i(o))^/^ (42) 

a known result, while the contact at unitarity, 

a strange looking result, but quite interesting, since this relationship and (142|) should 
be valid in general not only at the mean-field level. From the equations above one 
can check that these results follow by assuming universality at unitarity, LDA and the 
boundary condition p(r) — as r — )■ Vmax = Mean-field provides the 

values /i(0) and c(0) at unitarity, but these could be taken from experiments or more 
precise theoretical calculations. At mean-field level, we have shown that c(0) ~ A^(0), 
and therefore, knowledge of the contact in the trap at unitarity yields information on 
the uniform gap at unitarity. 

Given the approximated nature of mean-field, the calculated reduced contact in 
the trap c{fj) is in reasonable agreement with the more exact calculations by Werner 
et al.[3l] and with the recent experiments by Partridge et al.[9] and Stewart et al.[23]. 
This is seen in Fig. HI c vs ?), where we also plot experimental data from Ref. [9], as 
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taken from Ref.[34j, see their Fig. la. In Ref.jH] the number of bound pairs in the 
closed-channel was measured and in Ref. [31] this number was related to the contact; 
the experiment was performed in a gas of ^Li confined by a harmonic trap. On the 
other hand, in Ref. |23] a comparison of experimental contact data is made with the 
theory of Ref. |34] . showing good agreement as well. Thus, a similar agreement with 
mean-field results can be inferred by comparing them with the theoretical results of Fig. 
la of Ref. [31] . Certainly, at unitarity mean- field does not quite agree as expected, but 
the theory rings true as fas as the overall physics is concerned. 

The other interesting results are the asymptotic values in the inhomogeneous case. 
In the extreme molecular BEC these are, 

/ 5 \'/' 

fx ^ —7]"^ + — - , c ^ 2fi as fj +00. (44) 
\32t]J 

Restoring the units for /i, it yields to lowest order, /i ^ —h^/2ma^, again, half the 
binding energy of one molecule [28]. For the BCS side we find, 

/i ^ 1 , Inc ^ Txfj as 17 — 7- —00, (45) 

Comparing with the corresponding equations for the homogenous case, eq. (l3T]) and ( 132|) . 
we find that c, c, fi and jl have the same dependence on their corresponding reduced fj 
or 17 to leading order. For the chemical potential, the equality between the homogenous 
and the trapped case is expected, namely, it should be true that ji ^ 01 ji ^ Sp 
the BCS side, while /x —h?l2ma^ in the BEC extreme, regardless of the confining 
potential. However, the situation is somewhat unexpected for the contact. That is, 
one finds that asymptotically at BEC and at BCS, the reduced contact is related to 
the reduced inverse scattering length in the same functional form for both cases. In 
the uniform gas, the contact is directly related to the gap, see eq.(l271). however, in 
the inhomogeneous trapped case the gap should be a local function and, therefore, the 
contact in the trapped case is an integrated square gap, namely, 

C(V,/i,r/) = £^/A2(f)rfV (46) 

Even though the gap is a local quantity for a gas in a trap, indirect measurements of it 
have been alluded to in Refs.[9l [161 122], without fully considering its local nature, but 
certainly showing consistency with theory and their own data interpretation. Although 
this certainly needs a more careful analysis, we believe the reason behind this consistency 
may be both, the fact that one may ascribe a kind of "mean" gap in a trap, and that 
the dependence of the contact with r] is very similar in the homogenous and the trapped 
cases. 

To end this section, and because this suggests an independent and useful 
experimental determination, we calculate the corresponding generalized pressure for 
the trapped gas. As discussed at large in Refs. [lOl HH |12] , the grand potential in the 
trap is f2 = —W thus defining a generalized, pressure P, conjugate to the generalized 
volume V. This generalized pressure is the bona-fide counterpart of the hydrostatic 
pressure p of the homogenous gas. The equation of state of the trapped gas can be 
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Figure 4. Reduced contact c vs r) and chemical potential fi vs f), for a fluid in a 
harmonic trap, fj — I/kfo,- The data points are from the experiments of Ref.[9], as 
taken from Rcf . ^34^ . 

given in terms of V = V{N/V,T,ri) just as j9 = p{N/V,T,r]) in the homogenous one. 
It turns out that the measurement of V appears as a relatively simple task, given the 
enormous recent experimental advances. However, to make a direct comparison with 
the homogenous case, we will calculate the reduced grand potential a) instead. One can 
show that uj = —3V. To perform this evaluation we could resort to integration of the 
homogenous grand potential density using LDA. However, as shown in Refs. [10 | 1 ^ H2] . 
one can instead use an extremely simple formula for the grand potential in the trap. 
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Figure 5. Reduced generalized pressure p vs fj for a homogenous fluid, fj = l/npa. 
We have plotted 3p to make direct comparison with the homogenous pressure, see 
Fig El See text for details. 



namely, 



2 f ^1 

Q = — / p{r)-muj'^r'^ dPr. 

3 ^ 2 



(47) 



That is, sole knowledge of the density profile and the external potential suffice to know 
the grand potential. We insist, if one has simultaneous knowledge of the number of 
particles iV, the temperature T, and the scattering length a, then the equation of state 
of the system is immediately obtained. Using T = and reduce units for the trap, we 
find the following formula 

,5/2 



-647r 



X p{x; fj) dx, 



(48) 



, (A + v'f 

where again the upper value is used for fj < and the lower one for fj > 0. The result 
is shown in Fig. [5l Comparing with Fig. Ej we find once more a remarkable similarity 
with the homogenous case. 



5. Final Remarks 



In this article we have performed a simple and straightforward derivation of the contact 
C as the thermodynamic conjugate variable to the inverse of the scattering length a. We 
have shown that essentially all the known results regarding the contact follow from this 
identification. We see that this identification is independent of any microscopic model 
one uses, as long as the scattering length can be varied externally. For instance, the 



Thermodynamic origin of the Contact 



17 



state of the system may have any additional parametric dependence on other interatomic 
potential parameters and still the contact remains meaningful. Clearly, this is also valid 
for gases of Bose atoms. In the case where one can disregard other atomic parameters 
and keep the scattering length as the only two-body interaction dependence, dictated by 
the physical situation at hand, then a "universal" or law of corresponding states emerges; 
namely, in the appropriate reduced variables, all thermodynamic variables depend on 
the reduced temperature and inverse scattering length. 

Using the mean-field BEC-BCS theory, as introduced by Eagles [27] and Leggett[28]. 
one can explicitely exemplify all the above description. Furthermore, one finds a yet 
another identification of the contact, that is, its relationship to the (square) gap A. This 
may serve as an additional and independent way of measuring or determining qualities of 
the atomic excitations spectra. It is somewhat surprising that this approximated theory, 
in principle not valid in the crossover BEC-BCS region, still gives reasonable quantitative 
results, comparing well with both quantum Monte-Carlo calculations [5T| [52| \53\ ISU [55] 
as well as with actual experimental data[5l [HI [191 [201 |22l [23]. We would like to stress 
the fact that there is a remarkable similarity of the dependence of the reduced chemical 
potential, contact and pressure, on the reduced inverse of scattering length, for the 
homogenous and the trapped gases. This is not at all obvious from the outset and we 
wonder if this similarity is privative of the mean-field model or if it is a more general 
result. It is curious to note that, in this case of mean- field with sole dependence on 
the scattering length, the asymptotic values at the BEC and BCS extremes are more 
"universal" than at unitarity! It remains as a useful excersice to work out all the 
thermodynamics here presented using field theoretical extensions of mean field, see Refs. 
[56] and [57], to obtain more precise values at unitarity. 

In this article we have also pointed out the simplicity and, yet, the relevance of 
the measurement of the generalized pressure for trapped gases. Besides being readily 
measured just by knowing the density profile, the equation of state of the system is 
actually given in terms of this quantity, the frequency of the trap, the number of particles 
and the temperature. Its study and experimental determination may shed light on other 
properties, such as phase transitions and heat capacities [42j. 
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